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We study the effects of preparation of input states in a quantum tomography experiment. We 
show that maps arising from a quantum process tomography experiment (called process maps) 
differ from the well know dynamical maps. The difference between the two is due to the preparation 
procedure that is necessary for any quantum experiment. We study two preparation procedures, 
stochastic preparation and preparation by measurements. The stochastic preparation procedure 
yields process maps that are linear, while the preparations using von Neumann measurements lead 
to non-linear processes, and can only be consistently described by a bi-linear process map. A new 
process tomography recipe is derived for preparation by measurement for qubits. The difference 
between the two methods is analyzed in terms of a quantum process tomography experiment. A 
verification protocol is proposed to differentiate between linear processes and bi-linear processes. We 
also emphasize the preparation procedure will have a non-trivial effect for any quantum experiment 
in which the system of interest interacts with its environment. 
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I. INTRODUCTION 

Quantum process tomography [l|, [3| is the experimen- 
tal tool that determines the open evolution of a system 
that interacts with the surrounding environment. It is 
an important tool in the fields of quantum computation 
and quantum information that allows an experimenter to 
determine the action of a quantum process on the system 
of interest. 

The standard tomography procedure and some vari- 
ations of it (namely entanglement assisted tomography 
ancilla assisted tomography [^[^, and direct char- 

_) have been veri- 

II, [Hi. 




acterization of quantum dynamics 
fied experimentally [ol. [lol.]lll. \vX \Vc 

In many of these experiments, the maps that charac- 
terize the quantum process have been plagued with neg- 
ative eigenvalues and sometimes non-linear behavior. It 
was pointed out previously that the negative eigenvalues 
may be due to the initial correlations between the system 
and the environment [H, [13, [H, [H, [l^l . 

Quantum process tomography is thought of as a pro- 
cedure that allows us to experimentally determine the 
dynamical map describing the quantum process. How- 
ever, all experiments require a method to prepare the 
initial states of the system at the beginning of the exper- 
iment. This act of preparation has been neglected from 
the theory of quantum process tomography. 

We will show in this paper how an open system is 
prepared into different initial states can fundamentally 
change the outcome. A consequence of our study will 
force us to incorporate state preparation into the map 
that describes the process. We will have to distinguish 
between maps determined from a quantum process to- 



mography experiment, which we will call process maps, 
from the well know dynamical maps pll . [2^ . The key 
difference between dynamical maps and process maps is 
that process maps include the initial step of state prepa- 
ration, while dynamical maps are not restricted in that 
sense (see Appendix IB]) . 

We will study two methods for preparing states 
for quantum experiments, the stochastic preparation 
method and preparations using von Neumann measure- 
ments. In the former case, the process is given by a 
completely positive linear map (also see (20| for an inde- 
pendent, but related discussion). However for the mea- 
surement method, we will show that the outcome of the 
experiment cannot be consistently described by a linear 
map. We propose a bi-linear process map to describe 
such a quantum process, and we will show that this bi- 
linear map can be experimentally determined by devel- 
oping a procedure for bi-linear quantum process tomog- 
raphy. 

Based on our results existing quantum process tomog- 
raphy experiments should be reanalyzed. Any experi- 
ment which obtained a process map that had negative 
eigenvalues or behaved in non-linear fashion may suffer 
from poor preparation of input states, and should be an- 
alyzed for bi-linearities. 



II. LINEAR QUANTUM PROCESS 
TOMOGRAPHY 



The objective of quantum process tomography is to 
determine how a quantum process acts on different states 
of the system. In very basic terms, a quantum process 
takes different quantum input states to different output 
states: 
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The complete behavior of the process is known if the 
output state for any given input state can be predicted. 

The tomography aspect of quantum process tomogra- 
phy, is to use a finite number of carefully selected input 
states instead of all possible states, to determine the pro- 
cess. If the quantum process is described by a linear map 
(see [1^ for detailed discussion) , then the necessary input 
states should linearly span the state space of the system. 
For a finite dimensional state space, this requires a finite 
number of input states. Once the evolution of these in- 
put states is known, then by linearity the evolution of 
any input state is known. 

For example for a qubit only the following four projec- 
tions as input states are necessary 

= 1(1 + ^2), and P(3,+) = 1(1 + ^3) (1) 

to linearly span the whole state space, i.e. any state of 
a qubit can be written as a unique linear combination 
of these four projections. Above, 1 is the 2x2 identity 
matrix and aj are the Pauli spin matrices. 

Using the set linearly independent input states P^^\ 
and measuring the corresponding output states Q^-'\ the 
evolution of an arbitrary input state can be determined. 
Let the linear map describing the process be given by A, 
and the arbitrary input state be expressed (uniquely) as 
a linear combination PjP^-'^ ■ Then the action of the 
map in terms of the matrix elements is as follows: 



r's' \ j I j 

The map itself can be expressed as 

A , , - Vc>(")p(") * 

■"irr'.ss' — / , ^rs ^r's' ' 



(2) 



where P*^"^ arc the duals of the input states satisfying 
the scalar product 

p(m)tp(n) _ p(m)*p(n) _ r 

11— / ^ 1 1 ,.g — Ujyin- 

rs 

The duals for the projections in Eq. ^ are 

P(l--) = 1(1-^1-^2 -(73), 

P(l-+) = 1(1+^1-^2-^3), (3) 

p(2.+) = cr2, and P(3.+) = 0-3. 

III. PREPARATION OF INPUT STATES 

It is important to note the assumption made at the 
beginning of the last section, that is, the process is given 



by a linear map. This assumption is made without any 
consideration to how the input states are prepared. We 
will show that for the same input state, evolving through 
the same Hamiltonian, the output state would be dif- 
ferent depending on whether the input state is prepared 
stochastically or by a measurement. 

The basic steps in a quantum process tomography ex- 
periment arc broken down below: 

i. Just before the experiment begins, the system and 
environment is in an unknown state, which wc will 
write as 70. The system and environment could be 
entangled or correlated. For our discussions we will 
label the system as A and the environment as B. 

a. The system is prepared into a known input state. 
Let ^(") be the map that prepares the system into 
the n*^ input state. The system and environment 
state after preparation is therefore given by (^(")(g) 
X)(7o), where X is the identity map acting on the 
space of the environment. 

Hi. The system is then sent through an unknown quan- 
tum process. We consider the evolution to be a uni- 
tary transformation U in the space of the system 
and environment: 

f/(^(")(8X)(7o)C/t. 

iv. The trace with respect to the environment is taken 
to obtain the output state of the system 



Q 



(") 



Trm 



[/(^(") (^X)(7o)[/^ 



(4) 



where Q^"-* is the output state corresponding to the 
input state prepared by 

V. Finally using the input and the output states, a 
map describing the process is constructed. 

It is important to keep in mind that these basic steps 
(excluding the last step) also describe most quantum ex- 
periments, not just specifically quantum process tomog- 
raphy experiments. Therefore the following results may 
be applicable to many quantum experiments, not just to 
quantum process tomography experiments. 



IV. STOCHASTIC PREPARATION 

There are two methods in quantum theory to prepare 
an unknown state into a known state; stochastic prepara- 
tions and preparations by measurements. In this section, 
we will discuss the method of using stochastic maps for 
preparation of states. 

Let us consider preparation by a stochastic pin map 
which maps all states to a fixed single state i2J|. Using 
a pin map for preparation is a common initial step in 
various experiments. For example, in quantum dot ex- 
periments the system is cooled very close to absolute zero 
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temperature. This ensures the probabihty of the system 
being in the ground energy state is nearly one. This is 
effectively a pin map to the ground energy state. 

The experiment procedure will begin with a pin map 
Q, which takes any density matrix to a fixed pure state 
I $) . Then the state of the system and environment after 
the pin map is: 



(e®x) (70) = 1$) ($|®T(e), 



(5) 



where t{Q) = TrA[(8 (8) T) (70)]- The pin map fixes the 
system into a single pure state, which in turn means that 
the state of the environment is fixed into a single state 
as well. The purpose of the pin map is to decouple the 
system from the environment, to eliminate any correla- 
tion between the system state and the environment state. 
Note that the state of the environment does depend on 
the choice of the pin map 8. 

Once the pin map is applied, the system has to be 
prepared into the various different input states for the 
tomography experiment. This can be expressed in the 
most general way with stochastic maps: 

(1$) ($1) ^ (6) 

where P^"^ are the desired input states. 

The preparation procedure can be summarized as 
= fiM o Q. The overall experiment can be written 
in a single equation by combining Eqs. ([IHB]): 



Q(") = TrB 



e 



(70) t/t 



= Trn [C/P(") ®T(e)C/^ 



(7) 



We will call this equation 'the process equation'. No- 
tice that this process equation is linear on P^"^. Once 
the input states arc prepared, the procedure for quan- 
tum process tomography is the same as given in section 
im This is a generalization of linear quantum process 
tomography. 

It should be emphasized that the initial pin map O is 
critical; because for the process to be linear the state of 
the environment must be independent of the input state. 
It may be tempting to simply use a set of pin maps, 6^"\ 
to prepare the various input states P^"^. However, the 
process equation in this case yields: 



Q 



(n) 



Trm 



Trm 



c/e(")(g)X(7o)t/^ 
c/p(") ®T(e("))c/t 



(8) 



This is no longer a linear equation on p(") , since the state 
of the environment, r(0'"^), is effectively dependent on 



A. Example of Stochastic Preparation 

It may instructive to look at a simple example involv- 
ing two qubits at this point. We will treat one qubit as 



the system of interest and the other as the unknowable 
environment. 

Consider the preparation by a pin map O 



e®x(7o) = |$) ($1 



1 



-1. 



(9) 



that yields a pure state \^) for the system qubit and a 
completely mixed state for the environment qubit. 

The next step is to create different input states using 
different maps In this case, the fixed state 1$) ($| 

can simply be locally rotated to get the desired input 
state P("^ given in Eq. ([T|) 



p{n) 



:1, 



(10) 



where n = {(1, -), (1, +), (2, +), (3, +)} and are the 
unitary operators acting on the space of the system. Now 
each input state is sent through the quantum process. 
The output states can be calculated in a straight forward 
manner using the process equation Eq. ([7]). For this 
example we chose the unitary transformations U to be 



U = e 
The output states are 
1 



(11) 



1 



= i(l + CV2), and Q(3-+)^ 1(1 + ^^3). 



The linear process map is constructed using Eq. ([2]), 
the duals in Eq. ([3|), and the output states: 



A. 



(l + C 



y 2C^ 



2 










2C2 


1-C2 
1 + C2 



where C = cos(2i). 

We can verify that the process is correctly represented 
by this linear map by calculating As(/o) and comparing 
it to a direct calculation Tr^lUp® ^tU^, for any p. 



V. PREPARATION BY MEASUREMENTS 

Making a quantum measurement is another effective 
method of preparing input states for an experiment. For 
our discussion, we will use the von Neumann model 
for the measurement process. With the von Neumann 
model, a measurement is given by a set of orthonormal 
projections. If a particular outcome is observed from the 
measurement, the state of the system collapses to that 
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corresponding projection. Therefore, the input states can 
be prepared for an experiment by suitably fixing our mea- 
surement basis. With the knowledge of the basis and the 
outcomes, the exact input state is known. 

For experiments dealing with polarization states of 
photons this is a very practical method for preparation. 
Sending a randomly polarized beam of photons through a 
polarizer will yield in a perfectly polarized beam of light. 
This polarized beam can e used as one of the input state 
for a quantum tomography experiment. 

Suppose the n*'' input state is prepared, given by the 
projection by measurement. For the open system, 

the state becomes 

® X(7o) = ^ (P("' ® a) 70 [P^""^ ® l) , (12) 

where F^") — Tr[(p("^ l) 70] is the normalization fac- 
tor (it is the probability of obtaining that particular input 
state from a von Neumann measurement). For simplic- 
ity, we will from this point on, simply write p(") instead 
of P(") (g) 1. 

Combining Eq. ^ and Eq. (fT^ leads to the process 
equation that will relate the output states to the input 
states. For the n*^ input state the process equation is 



Q 



(n) 



-Trm 



[/P'"'7oP("^t/^ 



(13) 



Is this process given by a linear map? Dynamically, 
the evolution of the total state 70 is linear. However, for 
the purpose of tomography, the dynamics of the prepared 
input states P'^"^ is of interest, but P^"-* appears twice 
in the process equation, therefore the output states Q^") 
depend bi-lincarly on P'"'. That is suppose P — a A + 
(3B, where A and B are some operators and a and (3 are 
constants then 

P70P = (aA + /3P)7o(a^ + /3P) 

= a^A^oA + a(iA'^oB + a(iB"foA + l3'^B"foB. 

The last result shows that the Q has a bi-linear depen- 
dence on P. This bi-linearity can also be seen from the 
dependence of the environment state on the state of the 
system. To see this expanding Eq. (|13p as: 



= Tr™ 



[/p(") ^^(«)[/t 



with 



r(") 



TrA 



P(")70 



The last equation clearly shows that the environment 
state T^"-* is in effect a function of P^"-'. It is well known 
[19I . [20l . [2^ that if the initial state of the system is related 
to the state of the environment by some function /: 



P = 



where where is the state of the system and f{p^)* 
are the density matrix in the space of the environment. 



Then the evolution of the reduced matrices p cannot 
be consistently described by a single linear map. In this 
case, the function / is of a specific form that gives us a 
bi-linear dependence. 

We now look at a simple example to demonstrate that 
when input states are prepared by measurement, the re- 
sults cannot be consistently described by a linear map. 



A. Example of Preparation by Measurements 

Our example demonstrates when the output states are 
bi-linearly related to the input states the standard quan- 
tum process tomography procedure will give incorrect 
results. Suppose we are unaware of the bi-lincar depen- 
dence that arises from preparing input states by mea- 
surements, and assume that the process is given by a 
linear map. We would prepare a set of linearly indepen- 
dent input states, then construct the linear process map 
from the measured output states and the duals of input 
states. We will show that this linear map will not give 
the correct prediction for certain inputs. 

Consider the following two qubit state as the available 
state to the experimenter at i = 



70 



1 



a + C23Cr2 «i era) . (14) 



We will treat the first qubit as the system and the second 
qubit as the environment. Note that the the total state 
is separable, but it is correlated with the environment. 
Once again we will use the input states given in Eq. 
The state of the system plus the environment after 
each measurement takes the following form: 

^ pin) ^ 1 1 (fo, n^iil, ±), (3, +)}), 



p(2,+)^^p(2.+) ^ p(2,+) ^ 1 ( ^ 



C23 
1+02 



^3 



(15) 



Following the same recipe as before the output states 
are obtained using Eq. ([T^ . We will use the same uni- 
tary U as in the last example given in Eq. pT|) . The 
output states arc as follow: 



^1 



and 



- (1 - 4,CSa, + C'a2 + C2+35V3) , 
Q(''+) = ^(l + ^'^3), 



where C = cos(2<), S ~ sin(2t), and c^g = j+fj- 

The linear process map is constructed using Eq. ([2]), 
the duals in Eq. ([3]), and the output states: 



1 + 



■ C ZC23CS' 





2C^+ictCS 
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Now consider the action of on state 

p(2, 



o'2)- P^^' is a linear combination of the input states as 

p(2,-) ^ p(U-) ^ p(l. + ) _ p(2. + ) If tl^g j^c^jQj^ of jg 

hnear then the output state corresponding to the input 
state 

-) should be 



p(2 



)(2,+) 



(16) 



-23 



^23 



M fully describes the process, and we will call M the 
bi-linear process map. 

M contains both U and 70, however knowing M, is not 
sufficient to determine U and 70. As expected, it should 
not be possible to determine U and 70 through measure- 
ments and preparations on the system alone, without ac- 
cess to the environment. Conversely, M contains all the 
information that is necessary to fully determine the out- 
put state for any prepared input state. 

Before proceeding to determine the matrix M, let us 
look at some basic properties of M. 



Is this the same as if the global state was prepared in 
the state P^^^^) by a measurement? The output state for 
input 

p(2,-) 

can be calculated using Eq. p3p . 



Q 



(2,- 



1 



^(2,-1 



-Trm 



[/p(2.-)^„p(2.-)[/t 



- (a - c^,CSai - CV2 - 02-352^3) , (17) 



where Cnn = rr^^. 

Clearly the output state predicted by the linear pro- 
cess map in Eq. (jl6p is not the same as the real output 
state calculated dynamically in Eq. (|17p . hence the linear 
process map A^ does not describe the process correctly. 
This is not surprising, the state of the environment in Eq. 
(fTS]) depends on 02, and subsequently the linear process 
map depends on 02, hence on the state of the system. 
This is where the non-linearity of the process arises. 

In next section we will show that this process can be 
correctly described using a bi-linear process map. We will 
also develop a new tomography procedure for bi-linear 
process maps. 



VI. CHARACTERIZATION OF BI-LINEAR 
QUANTUM PROCESSES 

Let us expand Eq. with matrix indices: 



— TT p(") ^, p(") rr* 

ir,s ~ -p(n) '^re.r'a-'^r'r" Wr" a.s" (3-'^ s" s' ^ s<l,s' 



p(7 



^ Pj.iij.1 {Ur<L,r' alQr" a,s" se,s' 13) P< 



(n) 



p(n) ^ r"r' -^•-^j-" -^s" s' s" s' 

In the last equation, the matrix M is defined as: 



(18) 



Q/3e 

Note that in Eq. [TH]the superscript indices on A/ match 
the elements on the left hand side of the equation, while 
the subscript indices are summed on the right hand of 
the equation. 

The output state Q'") is given by the matrix M acting 
bi-linearly on the input state P'"^ Therefore the matrix 



A. Some Basic Properties of M 

Let us start with the trace of M. To take the trace 
equate indices r with s, r' with s' , and r" with s" to get: 

Tr[M] = Y: ^^-Vt'v 

rr'r" 

— Ure,r'a^0r"a,r"l3'^re,r'p- 
afJe rr'r" 

Since [/t[/ = I ^ Ureyp^reya = Sap, then 

Tr[A/] = J2^°r"a,r"a = 1- 

ar" 

As with the case of linear process maps, matrix M is 
hermitian. This is easy to see by taking the complex 
conjugate of M, 

^M^„'^,'^„^, j = j ^ t^re,r'Q70r"a,s"/3^s*£,s'/3 1 
\at3e J 

~ ^ Use,s'l3l0s"p,r"aUrt,r'a, 
apt 

The complex conjugate of M is not only the transpose of 
M but each element of M is also transposed. 



B. Bi-linear map for a qubit 

We now have to develop a new tomography procedure 
to deal with the bi-linear process map M . We will need 
to figure out a finite set of input states and the corre- 
sponding output states that will allow us to determine 
M . Afterall, this is the objective of tomography - by 
performing measurements on a small number of select 
input states a complete description of a process can be 
obtained and predict the output state for any arbitrary 
input state. 

M is a large [N^ x N^) matrix, where N is the dimen- 
sion of the space of system. To make it more manageable, 
Eq. can be interpreted as: 



(p(")|M|p(")) = r(")Q("^ 



(20) 
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where P^") is now treated as a vector. In this form, M is a 
N"^ X N"^ matrix with each element being a. NxN matrix. 
We will call Eq (|20|1 the ^bi-linear process equation\ 

Since M is a Hermitian matrix, it has ■^{N'^ + N"^) 
independent block elements. Therefore ^{N'^ + N^) in- 
dependent equations in the form of Eq. (PP]) are neces- 
sary to fully determine M. It is clear that neither an 
orthonormal set of N or linearly independent set of N'^ 
input states would provide sufficient number of equations 
to resolve the bi-linear process map M. 

Consider a qubit system. In that case the projections 
P'"' can be written in terms its coherence vector aj and 
Pauli spin matrices a^: 



p(r. 



E 



(") 



(21) 



Since P*^") is a projection then there is an additional 
constrain = 1- 

The matrices 1 and Uj together forms a vector basis 
for this space. Therefore Eq. pT|) is simply a vector 
decomposition of P^") in a fixed basis. Taking this form 
for p(") and substituting into Eq. (pn|) gives: 



r(")Q(«) = (p(")|A/|p 



(IIMID 



(") (") 



J") 
J— 
4 



((l|M|a,)-f (a,|A/|U)) 



a,|A'/|a,) 



(22) 



Observe that the terms (1|M|1), {t\M\aj), (crj|A/|l) 
and {aj\M\ak) , are simply the matrix elements of M in 
basis. We just need to find a set of projections 
P'-"-' that will allow us to solve for these matrix elements. 

Consider the following specific projections defined as 
pU,±) = i(l±crj) with.? = {1,2,3}. 



p(j,±)gb-,±) ^ ^pO-,±)|^,^|pb-,±) 



-mM\l) + {a,\M\n,)) 
±i((a,|M|a) + (a|M|a,)). 



(23) 



Simultaneously solving the {+) and (— ) equations above 
give the following unknowns: 



To obtain the cross terms {a.j\M\iJk) consider projec- 



tions such as p(-'+ 



fc+i. 



1 



T^c^k) for 



k > j which gives: 

p(j+fc+l,+)g(i+fe+l. + ) ^ |p(i+fc+l- + )^ 

= i((l|A/|a) + (a,|A/|a,)) 
+ i((l|A/|l) + (a,|A/|f7fe)) 
^ ((l|A/|a,) + (f7,|A/|a)) 

((l|Af|(7fe) + (f7fe|M|a)) 



4^/2 
1 



+ -((a,|Af|afe) + (a,|Af|a,)). 

Substitute the known terms and solve for the desired 
cross terms, 

{a,\M\ak) + {<7k\M\a,) = -2 (l + ^/2) r(^"-+)g(^"-+) 

-2(l-\/2)r(J'-)Q(^^-) 

-2(i + V2)r('^-'+)g('=^+) 
-2 (i - V2) rC^-'-^QC^---) 

In summary using the following nine projections, 

and solving them simultaneously yields all desired matrix 
elements: + {aj\M\aj), {t\M\a.j) + (crj|Af|l), 

and {(Jj\M\ak) + {ak\M\aj) . 

Though this is not enough to fully determine M, these 
elements are sufficient to determine the output state for 
any input state. Using the property ^ 1' Eq. 

can be rewritten as: 



+ 5]a5")((l|A./|a,) + (a,|A/|l)) (25) 



+ J2''j''k(.{'^j\M\a,) + {a,\M\a,)). 

Observe that the sums of the cross terms {aj\M\(Tk) + 
{ak\M\aj) can appear together because the coefficients 

Cj"^ are real. Also, the element (a|Af|l) can always be 
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paired with a diagonal clement {aj\M\aj) as long as the 
state is a pure projection satisfying = 1- The 

diagonal element (l|Af |1) only has to be known if the 
system can be prepared directly to a mixed state such 



that 



,(")\2 



< 1. This may be accomplished by a 



generalized measurement [2g| (see appendix 1^ . If gen- 
eralized measurements are allowed, then just one more 
input state is needed, for example ^(1 + ^cti), which 
gives another independent equation that can be solved 
to obtain (1|M|1). 

Therefore the elements of M found in Eq. (pS)) are 
all that are needed to describe the process. By measur- 
ing the outputs for the nine specified input states, the 
matrix M can be calculated. We now have a good quan- 
tum process tomography procedure for an open two-level 
system. 

This procedure can be generalized for an n-lcvel system 
in the same spirit as above. However care must be taken, 
since the constraints on the coherence vector of an n- 
level system are more complicated (see [13] for a detailed 
discussion) . 



VII. BI-LINEAR VS. LINEAR PROCESS MAP 
VERIFICATION PROCEDURE 

In practice it is not easy to tell if a process is linear 
or bi-linear. The bi-linear process map is incompatible 
with the behavior of a linear process map and conversely 
a linear process map will not adequately describe a bi- 
linear process. We have developed a procedure below 
that will differentiate between a linear process and a bi- 
linear process. 

Consider what happens to a state that is a linear com- 
bination of a set of projections: 



If the evolution of X though the process is bi-linear then 
the output is written in terms of the bi-linear process 
map as: 



Therefore the bi-linear process map gives different predic- 
tions from a linear process map. In that case, we should 
be able to distinguish between whether a process is given 
by a linear map or a bi-linear map. 

Consider the example of the tomography procedure 
proposed for a two-level system. If the process is given 
by a linear process map, then the nine input states in Eq. 
(j24p are over complete; only four input states are needed 
to determine a linear process map. This discrepancy is 
summarized by the following linear sum rules: 

p(l ,-) _ p{2.+)^ 
p(l -) _ p{3,+) 



x/2 



p(2,-) 




p(3,-) 


_ p(l,+ 


p(4,±) 




p(5,±) 


^^(- 


p(6,±) 


-hi"' 



± 1 fp{3,+)_p(l,-) 



(1,+) +p(i' 

(2, + ) ^ p(3, + ) _ p(l, + ) _ p{l,- 



If the process is linear, then the output states must 
satisfy the same sum rules, which are obtained from the 
above equations by suitably writing Q in place of P. If 
these sum rules are not satisfied, then the process is not 
linear. 

However satisfying the sum rules is necessary but not 
sufficient to determine if the process is linear; a bi-lincar 
process map can still be constructed from this set of input 
and output states without contradiction. Therefore an 
additional input state, distinct from the above nine input 
states, should be tested to determine which of Eq. ([77|) 
and Eq. (|28p is satisfied. 

More explicitly, if the output state corresponding to 
an additional input state of the form of Eq. 126]) is 
(26) found to be E„ c^r(")Q("), then the process is linear. 
However if the corresponding output state is given by 
^mn'^rnCn (P^'"' | A/ ) tlicn the proccss is bi-lincar. 



(X|M|X) ^ c™c„ (P^") I A'/ IP^-^ 



VIII. FUNDAMENTAL ISSUES OF STATE 
PREPARATION BEYOND QUANTUM PROCESS 
(27) TOMOGRAPHY 



But if the bi-linear process map is compatible in some 
way with a linear process, then Eq. (P7|) can be simplified 
to 



{X\M\X)^Y^ cl (P(") I M |P(" 



(28) 



In that case the two equations above are equal: 

Y^CmCn (P^"'! A/ |P(")) = (p(") |A/|p(")) . (29) 

mn n 

It is clear that no non-trivial conditions exists for M 
that will allow this equality for arbitrary coefficients Cm- 



Although our discussion has largely focused on quan- 
tum process tomography, we would like to emphasize that 
the process equations describe any quantum experiment. 
Before any experiment begins, the quantum system or 
particle would exist in an unknown state that could be 
(and most likely is) correlated with a quantum environ- 
ment. Preparation of the system or particle into a known 
state is a necessary part of any experiment. 

As we have shown, when the system or particle is 
imperfectly isolated, in other words, it interacts non- 
trivially with a quantum environment during the exper- 
iment, the initial step of state preparation in an experi- 
ment is of fundamental importance. 
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We described two methods of state preparation, the 
stochastic preparation and the preparation by measure- 
ments. For preparation by measurements, the outcomes 
of the experiment will be non-linearly related to the 
prepared states. So it would seem that the stochastic 
method is preferable, since with the stochastic method, 
the evolution of the prepared states to the final states is 
linear. However, how the stochastic maps arc actually 
performed may need to be carefully considered. 

Stochastic maps can be equivalently performed by a 
unitary transformation with the addition of ancillary 
quantum systems. In particular, the way the stochas- 
tic preparation is constructed in section [IV] is equivalent 
to a generalized measurement (see appendix E)) . 

Recall that with the stochastic preparation, different 
maps prepare different input states. 

Let K be the number of input states. Assume that for 
the experiment, an equal number of each input state is 
prepared, so the expectation map is: 

K 



n=l 

This overall, expectation map is a trace preserving map 
since each preserves trace. Therefore this can be 

considered to be a generalized measurement, where each 
^^{n) represents a measurement outcome. 

A generalized measurement can be equivalently per- 
formed by a unitary transformation and a von Neumann 
measurement by using a suitable ancillary system. How 
docs this implementation of the stochastic preparation 
effect our results? Let the ancillary system be labeled as 
C, and the generalized measurement is implemented with 
a unitary transformation W and von Neumann measure- 
ment given by the projections J*^"^; where the J^"^ out- 
come corresponds with the n^^ preparation map 



B, and the ancillary C then we cannot make this simpli- 
fication, and Q^"^ would have a bi-linear dependence on 
, and in turn a bi-linear dependence on the prepared 
input states. 

We can make two important observations from this. 
First, the quantum environment should be defined as any 
quantum system that interacts with the primary system 
A during the experiment. The quantum environment 
does not have to include any quantum systems that the 
primary system is entangled or correlated with, if there 
is no interaction between the two during the experiment. 

The second observation is that if the stochastic prepa- 
ration method is used, any ancillary systems used to im- 
plement the stochastic maps, must be perfectly isolated 
from the primary system during the experiment. In other 
words, they must not be part of the quantum environ- 
ment. If the ancillary systems used are not properly iso- 
lated, then the process may have bi-linear dependence on 
the input states. 

This is an important result - simply deciding on a 
stochastic preparation method for an experiment may 
not in practice guarantee that the process will be linear. 
This can have implications in many areas, in particular, 
quantum error correction protocols are designed to cor- 
rect for linear noise, correcting for bi-linear effects may 
be a more difficult challenge. 

The verification procedure discussed in the last sec- 
tion may therefore be important, since it can be used as 
a tool to confirm the proper isolation of the apparatus 
and ancillary systems during the experiment. We can 
verify in practice if a process is linear instead of simply 
making that assumption. In the following section, we will 
describe a complete experiment that may be performed 
with qubits, including verification steps to check if the 
process is linear or bi-linear. 



(30) 



where is the initial state of the ancillary system. 

Now include the above equation within the overall pro- 
cess equation Eq. ^ 



IX. A COMPLETE RECIPE FOR AN 
EXPERIMENT 



If U acts only on the system A and environment B, 
and not on the ancillary C, then the above equation can 
be simplified: 



l(n) 



Tri 
Tri 
Tru 



t/Trc 



(31) 



And the result is the linear process map as before. How- 
ever, if U acts non-trivially on our system A, environment 



Although we have established the theory, let us make 
the ideas more concrete by developing a complete recipe 
for an experiment that can be used to determine whether 
a process is bi-linear or linear. We will also show specifi- 
cally how the corresponding bi-linear map or linear map 
can be calculated from the measurement results. 

For bi-linear process tomography, nine input states are 
necessary. For the nine states derived in section FVlB I the 
first six states are three pairs of orthonormal projections, 
but the last three are not. Now consider twelve projec- 
tions, nine from Eq. (|24p and three orthogonal to the 
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last three in that equation: 



p(4,- 
p(5,- 
p(6, 



1 - 



1 - 



1 


1 \ 























(32) 



These twelve projections are neatly grouped into six 
different sets of orthonormal pairs. If the states are pre- 
pared using von Neumann measurements, these would 
correspond to measurements in the ui, (72, cs, ci + (72, 
CTi + (73 and (72 + CTs directions. Twelve states are more 
than necessary for bi-linear process tomography, but the 
extra states can utilize as consistency checks. 

After recording the corresponding output states for all 
twelve input states, the linearity of the process can be 
verified. If the process is linear then only four linearly 
independent input states are necessary to determine it. 
The other eight input states can be written in linear com- 
bination of these four. If the states from Eq. ^ are used 
as the four linearly independent states then the following 
eight linear sum rules (that the input states satisfy) have 
to be satisfied: 



V2 



Q 



(6,±) _ 



^ V2 



)(2, + ) 



(3,+)_q(1, + )_q(1,-) 



If the eight sum rules are satisfied, then the process is 
not bi-linear, and should be described by a linear pro- 
cess map. The linear process map can then be computed 
following the recipe in section HIl 

If the eight sum rules are not satisfied, then the process 
may be given by a bi-linear map. We will attempt to 
verify that the process is bi-linear and calculate the bi- 
linear process map. 

Note that the probabilities T^") = Tr[7oP("'] associ- 
ated with each preparation should be found experimen- 
tally. The probabilities should be complete for an or- 
thonormal set of projections, in other words: T^^'~^^ + 
T^i-^^ = 1. Therefore the probabilities can be calculated 
from the fraction of the (4-) states as compared to the (— ) 
states, for all preparations made in the same direction. 

If the process is bi-linear then three additional states 
should evolve in a way that is consistent with a bi-linear 
map derived from the other nine states. The following 
equations are derived from this condition and using Eq. 
(PS)) . If these equations are satisfied then the process is 



bi-linear: 

r(4,-)g(4,-) 



(5,-) _ 



r(6,-)g(6,. 



^p(2,-)g(2,-)_p(2, + )g(2, + )^ 
^^(4,+)q(4. + )^ 

+r(3,-)Q(3,-)_r(3.+)Q(3.+)j 

1 /'p(2,-)g(2.-)_p(2. + )Q(2.+) 

v2 ^ 

^r(3,-)g(3,-)_p(3,+)g(3,+)^ 



If the conditions above arc satisfied then the process is 
bi-linear, and the bi-linear process map can be computed 
by following the recipe in section fVI Bl 

Once the matrix elements of M are determined in this 
fashion, the evolution of any state X = ^(1-1- J2jPj'^j) 
is given by: 

4rg = 4:{x\M\x) 

= J2p'^ {{t\M\l) + {<j,\M\a,)) 
j 

+ ^p,((l|A/|a,) + (a,|A/|l)) 

3 

k>j 

Note that since Q is a normalized state, the nor- 
malization constant T is the measurement probability 
r = Tr[X7o]. Although we had not explicitly mentioned 
this before, the matrix M contains all information about 
the measurement probabilities, that is why we needed 
the measurement probabilities F^"' to calculate the ma- 
trix M. 

Finally note that if both the test for linearity and bi- 
linearity fails, then the process cannot be consistently 
described by either a linear map or bi-linear map. The 
experiment then should be carefully analyzed for prob- 
lems such as any non-linear dependence that may have 
been introduced if the input states are not accurately 
prepared, or if there is some dependence of 70 on the 
prepared state. 



X. ANALYSIS OF A QUANTUM PROCESS 
TOMOGRAPHY EXPERIMENT 



In this section we will analyze a quantum process to- 
mography experiment performed by M. Howard et al. 
[31 ■ Our critique will emphasize the importance of hav- 
ing a consistent theory of state preparation. 
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In this experiment, the system that is studied is an 
electron configuration formed in a nitrogen vacancy de- 
fect in a diamond lattice. The quantum state of the sys- 
tem is given by a spin triplet (S=l). Again we will write 
the initial state of the system and environment as 70. 

The system is prepared by optical pumping, which re- 
sults in a strong spin polarization. The state of the sys- 
tem is said to have 70% chance of being in a pure state 
10). Or more mathematically, the probability of is 
Tr[|<^) (01 70] = 0.7. 

Since the population probability is high, an assumption 
was made that the state of the system can be simply ap- 
proximated as a pure state From this initial state, 
different input states can be prepared by suitably apply- 
ing microwave pulses resonant with the transition levels. 
After preparation, the system is allowed to evolve, and 
the final states (density matrices) are measured. With 
the knowledge of the initial state and the measured final 
states, the linear process map that should describe this 
process is determined. 

It was found that the liner process map has negative 
eigenvalues, so the map was "corrected" using a least 
squares fit between the experimentally determined map 
and a theoretical map based on Hermitian parametriza- 
tion psj . while enforcing complete positivity. 

However, if we do not regard the negative eigenvalues 
of the map as aberrations, then we should consider the 
assumptions about the preparation of the system more 
carefully. The assumption about the initial state of the 
system is: 



70 



(33) 



This is in effect a pin map. Together with the stochastic 
transformation of the initial state into the various in- 
put states, this is identical to the stochastic preparation 
method discussed in section HVl 

It is clear that the pure initial state assumption is un- 
reasonable, given our knowledge now of how the process 
is sensitive to the initial correlations between the system 
and the environment. In effect the action of the pin map 
in this experiment is not perfect, and therefore the pin 
map can be ignored. Then the process equation is: 



(34) 



where fi*^") is the stochastic mapping corresponding to 
preparing the n*'Hnput state. We will assume that the 
stochastic process does not involve any ancillary systems 
that interact with the primary system during the exper- 
iment. 

In this experiment, O^"-* is nothing more than a unitary 
transformation F^"^ satisfying V^^'^ = l^/'^"-*), where 
l^/i^")) is the desired pure n^^ input state to the process. 

We can write the unitary transformation for a two-level 
system as: 



where (^'^"^IV'i"^) = and (0|0±) = 0. This basically 
defines V^"^ as a transformation from the basis {\4>)} to 
the basis {l"*/"!"^)}- The equation for the process becomes: 



)(«) 



(35) 



TrB [t/|V^"^)(0|7o|'^) (V-^'^^lt/^ 
+TrB [t/|V'i"^)(0x|7o|0) (V^")|C/^ 
+TrB [c/|V^"))(0|7o|0i)(^^''|(7t 

+TrB t/|V'^^')(0x|7o|0±)(V'f 



Therefore, since (0|7o|0) = 0.7, to first approximation 
the process is a linear mapping on the states l^-*")) (-0'"^|- 
However it is clear that if all terms are included, the pro- 
cess is not truly a linear map of the states |?/;(")) (-0^"^|. 
The negative eigenvalues are therefore a result of fitting 
results into a linear process map when the process is not 
truly linear. 



XI. CONCLUSIONS 

Preparation of the system or particle into a known 
state is a necessary part of any experiment. We have 
shown that with open systems, some care has to be taken 
to define the method of state preparation. 

We described two methods, the stochastic method and 
the measurement method. The stochastic method leads 
linear processes that can be described by linear process 
maps. The advantage of the measurement method is that 
the primary system does not have to be perfectly isolated. 
This method only requires a good measurement appara- 
tus. 

With the stochastic method, the initial state can be 
made simply separable, effectively de-coupling the sys- 
tem from the environment. The evolution of the system 
is then given by a linear process map. However we find 
that the isolation of the apparatus from the system dur- 
ing the experiment is of greater importance with this 
method. Any apparatus or ancillary systems used for 
the stochastic preparation must not be contained in the 
quantum environment, the quantum environment being 
defined as everything that interacts with the quantum 
system during the experiment. 

The stochastic method is more consistent with the 
traditional method of performing quantum experiments, 
and its advantage is that the process map is effectively 
equivalent to the linear dynamical map. However the 
disadvantage to the stochastic method is that the appa- 
ratus and any ancillary systems employed to perform the 
stochastic transformations or generalized measurements, 
must be perfectly isolated from the primary system for 
the duration of the experiment. 

With the measurement method, if the initial state is 
correlated with the environment then the evolution of 
such a system is given by a bi-linear process map. The 
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determination of this bi-linear process map by process to- 
mography is more difficult, but we developed a procedure 
that works for qubit systems. 

If interaction occurs between the system and this envi- 
ronment, then non- linear noise can be introduced into the 
experiment. This could have fundamental consequences, 
for example quantum error correction schemes proposed 
so far are based on correcting linear noise. If bi-linear 
errors occur, error correction becomes a more difficult 
challenge. 

We proposed a protocol to distinguish between a bi- 
linear process and linear process. This protocol can be 
used to verify the assumptions made about state prepa- 
ration in the experiment, by adding some additional in- 
puts to the experiment and making consistency checks, 
the correctness of a stochastic preparation can be verified. 
The protocol we proposed can be an practical experimen- 
tal tool. 



The motivation is to use the most general linear op- 
eration on a density matrix as a measurement, a linear 
map. A general measurement can be framed as follows: 
given a quantum system in state p, let the measurement 
be given by a set of positive trace reducing linear maps 
B'-'^K The probability of registering the j*'' outcome is 
given by Tr[S'^-'(p)]. If the j*'' outcome is observed, then 
the system collapses to the state: 



Tr[SO-)(p)]- 



(Al) 



To ensure the probabilities sum to one. the sum of the 
maps yS^^' must be itself a trace preserving map. 

Writing B^^^ its canonical form with its eigen-matrices 
Ca^ and eigenvalues Cq ■* yields: 
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Appendix A: GENERALIZED MEASUREMENTS 



The condition that the overall map 'J2j ^^^^ preserves 
trace gives: 



(A3) 



Given this condition, we can prove that this very gen- 
eral scheme can be accomplished with a unitary trans- 
formation and von Neumann measurement, by using an 
ancillary system of sufficient size. The method we use 
here is similar to the method used to show how a trace 
preserving map can be given as a reduced unitary trans- 
formation 1321 . 



Overview of measurements 



Measurements in quantum theory begin with the von 
Neumann measurement p9j . which can be quickly sum- 
marized as follows: if the system being measured is in 
the state p and the measurement is given by a set of 
orthonormal projections 11^, then the probability of ob- 
taining the outcome is Tr[pnj]. If the outcome is 
observed, then the state collapses to 11^ . 

An early generalization of the von Neumann measure- 
ment was given [30,[3l| , with the introduction of POVMs. 
With POVMs, the measurement is still given by a set 
of operators Sj, however the operators need not be or- 
thonormal projections, but arc in general positive oper- 
ators satisfying "Y^^ = 1. The probability of obtaining 
the j*'* outcome is still Tr[pSj] and the state collapses to 

However, measurements can be generalized beyond 
POVMs (see [1^ for more discussions). Suppose a very 
elaborate measurement apparatus is constructed, one 
that involves many transformations and POVMs. Rather 
than to look at all the intricate details of this apparatus 
the apparatus can be treated as a black box. When a 
state is fed into the box one of several possible outcomes 
is registered, and a new state p' leaves the box, which is 
fixed for that particular outcome. 



2. Generalized measurements 

Consider a transformation W that acts on a tensor 
product space consisting of the original system and two 
ancillary systems. Let the original system and the an- 
cillary systems be spanned by the basis states \r,j,a) 
respectively. The action of W is defined as follows: 

W : |r', 0, 0) ^ E \f^\cli\r' |r, j, a) . (A4) 

rja 

The size of the ancillary systems is bounded by jiN"^ , 
since j ranges from 1 to p,, where p is the number of maps 
S'^-* making up the measurement and a ranges from 1 to 
N"^ (where A'' is dimension of the system) since each map 
S'^-* has at most N'^ C-matrices. 

The action of the transformation on a complete set 
of basis states is not yet defined. However for the states 
on which it is defined, W does preserve orthornormality 
between those states. The proof is straight forward using 
Eq. \M 

((r',0,0|W^t)(^|,'^0,0)) = E^"^[^"^]*v[Ci''^W 

rjoL 
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Since W preserves the orthonormahty, it can be made 
into a vahd unitary transformation by carefully defining 
its action on the remaining space that is so far not covered 
by equation IA4I 

Now we demonstrate that our generalized measure- 
ment, given by the set of maps B^^\ can be equivalently 
performed by this unitary transformation W and a von 
Neumann measurement. Performing the unitary trans- 
formation W on the original system in state p and the 
ancillary systems in the initial state |0,0) (0, 0| to gives: 

X = W{p(g)\0,0) (0,0\)W^ 

= E ^f^lC^'Kr'Pr'AC^p^rs's (A5) 

rsr' s' jka/S 

X \rj,a) {s,k,(3\ . 

Now perform a von Neumann measurement on the first 
ancillary system, given by the set of orthonormal projec- 
tions \j) The probability of the j*'' outcome is: 

Tr[x|j>(j|] = E ^^\C^^]rr'Pr's'\C^^]:'r (A6) 



Dynamical maps used to describe the evolution of open 
systems are usually defined with a constant environment 
state r: 

^(p*) =TrB [Up(^TU^ 
Implicitly, the state of the environment t is a parameter 
of the map Therefore, the linear dynamical map £i§ 
would only consistently describe an experiment if differ- 
ent input states p can be prepared independently of the 
environment state r. The actual issue of how this can be 
executed is never addressed. 

The stochastic preparation method provides a way for 
an experiment to be made so that different input states 
can be prepared with a fixed environment state. Consider 
the process equation Eq. ([7]) in section HVl 



g^") = Tn 

= TriB 



u 



o e 



= Tr[S(j)(p)]. 

If the j*'' outcome is observed, then the original plus 
ancillary system collapses to the state: 

ij\x\j) = i E ^¥^[Ci%r'Pr'AC^^^]:,s 



The process map is then given by: 



rsr' s'ap 

X \r,j,a) {s,j,P\ . 

where K = Tr[B^^\p)] normalizes the state. 

Finally the measurement is over and the particle exits 
the apparatus. The state of the system, now outside 
the apparatus, is given by tracing over the remaining 
ancillary system inside the apparatus: 



rs / J a 



l^a \ss' 



rsr' s' a. 
^ [B^-Hp)]rs 

Tr[6(")(p)]' 

Therefore this gives the same results as the general- 
ized measurement we laid out. We only needed an an- 
cillary system big enough (dimension pN'^), one unitary 
transformation, and one von Neumann measurement to 
perform the most general quantum measurement. 



Appendix B: PROCESS VS. DYNAMICAL MAPS 

Dynamical maps [U, [l^] allow the description of 
stochastic processes and the evolution of open systems. 



A(p-) 



Trm 



Up^®T{<d)U^ 



Therefore, in this context, the dynamical map is equiv- 
alent to the process map. However, for consistency, we 
have to remember that the environment state is a con- 
stant to the problem, therefore the pin map Q should 
also be a constant to the problem. 

It is possible to consider a dynamical map where the 
environment is not fixed, such as the reduced dynamical 
evolution of a non-simply separable state 70 [16| : 



[7o])=TrB[C/7ot^t] 



(Bl) 



The dynamical map in this problem is applicable only 
over a compatibility domain of states, rather than over 
the complete state space of the system A. The compat- 
ibility domain [l^, [33| is the set of states that are com- 
patible with the correlations in 70. Formally, this prob- 
lem defines an extension map [ssl. [3^. [35j (also known as 
preparations) that relates the initial state of the system 
A to the overall initial state of AB. The extension map 
is linear but not necessarily a completely positive map. 
Therefore, in this context, such dynamical maps are a 
theoretical tool, and have a limited relation to a physical 
experiment. This why we differentiate between process 
maps and dynamical maps. 
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